Vacuum friction in rotating particles 



A. Manjavacas and F. J. Garcia de Abajcj2] 
Instituto de Optica - CSIC, Serrano 121, 28006 Madrid, Spain 
(Dated: September 22, 2010) 

We study the frictional torque acting on particles rotating in empty space. At zero temperature, 
vacuum friction transforms mechanical energy into light emission and produces particle heating. 
However, particle cooling relative to the environment occurs at finite temperatures and low rotation 
velocities. Radiation emission is boosted and its spectrum significantly departed from a hot-body 
emission profile as the velocity increases. Stopping times ranging from hours to billions of years are 
predicted for materials, particle sizes, and temperatures accessible to experiment. Implications for 
the behavior of cosmic dust are discussed. 
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I. INTRODUCTION 

The radiation emitted by accelerated charges produces 
reaction forces acting back on them. For rotating charged 
particles (e.g., electric^ and magnetic^ dipoles) , this gives 
rise to reaction torquesP Likewise, accelerated neutral 
bodies are known to experience friction because they emit 
light due to the absolute change in the boundary condi- 
tions of the electromagnetic field. This is the so-called 
Casimir radiationPEl 

A spinning sphere presents a more challenging situa- 
tion: its surface appears to be unchanged, although it 
experiences a centripetal acceleration. So, the question 
arises, does a homogeneous, neutral sphere emit light 
simply by rotating? Is such a particle slowing down 
when spinning in vacuum? We know the inverse pro- 
cess to be true: the angular momentum carried by light 
can be transformed into mechanical rotation of neutral 
particlesP However, this type of problem requires a deli- 
cate analysis, somehow related to the non-contact friction 
predicted to occur between planar homogeneous surfaces 
set in relative uniform motion^ which is currently gen- 
erating a heated debateP 

In this paper, we investigate the friction produced on 
rotating neutral particles by interaction with the vacuum 
electromagnetic fields. Friction is negligible in dielectric 
particles possessing large optical gap compared to the ro- 
tation and thermal-radiation frequencies. For other ma- 
terials (e.g., metals), in contrast to previous predictions,^ 
we find nonzero stopping even at zero temperature. The 
dissipated energy is transformed into radiation emission 
and thermal heating of the particle, although cooling rel- 
ative to the surrounding vacuum is shown to take place 
under very common conditions. We formulate a theory 
that describes these phenomena and allows us to predict 
experimentally measurable effects. 

II. THEORETICAL DESCRIPTION 

We consider an isotropic particle at temperature T\ 
spinning with frequency f2 and embedded in a vacuum 
at temperature T (see Fig. [I]). The particle experiences 




FIG. 1: Sketch of a spherical rotating particle and parameters 
considered in this work. The particle is at temperature T\ and 
rotates with frequency Q. The interaction with vacuum at 
temperature To produces a frictional torque M and a radiated 
power P rad . 



a torque M by interaction with the surrounding radia- 
tion field and it is also capable of exchanging photons, 
with net emission power p rad . For simplicity, we assume 
the particle radius a to be small compared to the wave- 
length of the involved photons, so that we can describe 
it through its frequency-dependent polarizability a{uS). 
Since the maximum frequency of exchanged photons is 
controlled by the rotation frequency and the thermal 
baths at temperatures To and T\ , this approximation im- 
plies that both fla/c and kgTja/ch are taken to be small 
compared to unity. These conditions are fulfilled in very 
common situations (for instance, for a = 50 nm, one has 
fl < 6 x 10 3 THz and Tj < 4.6 x 10 4 K). 

Friction originates in fluctuations of both (i) the vac- 
uum electromagnetic field E fl and (ii) the particle po- 
larization p fl . We calculate the emitted power from the 
work exerted by the particle dipole, 

prad = _ ^ E i„d . Q p fl/ dt + E A . d P ind /dt) , (1) 

where E md is the field induced by p fl , and p lnd is the 
dipole induced by E fl . Likewise, the torque is obtained 
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from the action of the field on the dipole, 



M = (p fl x E 



ind 



,ind 



E fl ) 



(2) 



The result is quadratic in E fl for contribution (i) and in 
p fl for contribution (ii). The brackets () represent the 
average over these quadratic fluctuation terms, which we 
perform using the fluctuation-dissipation theorem (FDT) 
(see Appendix). 

Rotational motion enters here through the transfor- 
mation of the field and the polarization back and forth 
between rotating and lab frames. This is needed be- 
cause the particle polarizability can only be applied in 
the rotating frame, in which the electronic and vibra- 
tional excitations participating in a are well defined and 
il-independent. In contrast, the effective polarizability 
in the lab frame has a dependence on 57. Further details 
of this formalism are given in the Appendix. The result- 
ing radiated power reads (see Appendix for a detailed 
derivation) 
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is the spectral distribution of the rate of emission 
(when loT > 0) or absorption (ujT < 0), rij(uj) — 
[exp(Tiuj/kBTj) — is the Bose-Einstein distribution 
function at temperature Tj, 



gi(w) = Im{a;(w)} 



2uj 3 
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are odd functions of u> describing particle absorption 
for polarization either parallel (I —\\) or perpendicular 
(l =_L) with respect to the rotation axis, and p° = 
uj 2 /tt 2 c 3 is the free-space local density of photonic states. 
These results apply to particles with orthogonal princi- 
pal axes of polarization, rotating around one of them, 
and with a± given by the average of the polarizability 
over the remaining two orthogonal axes. The torque M 
takes a similar form, 



M = - 



duj hT(u>). 



(5) 



Incidentally, the term of Eq. Q vanishes under the 
integral of Eq. |5]), and furthermore, M = for Q — 
0. In the To = T\ = limit, one has rij(uj) = —0(u>), 
from which we find the integrals to be restricted to the 
(0, fi) range: only photons of frequency below f2 can be 
generated. 

Unfortunately, Eqs. ([I]) and ^ do not account for ra- 
diative corrections coming from the elaborate motion of 
induced charges in the rotating particle. Although such 



corrections are insignificant for small particles, we incor- 
porate them here for spheres in a phenomenological way 
through the term proportional to | a | 2 in Eq. Q, preceded 
by a coefficient chosen to yield gi = (and consequently, 
M = 0) in non-absorbing particles: 10 internal excitations 
(i.e., absorption) are necessary to mediate the coupling 
between the rotational state and radiation.EQ Further- 
more, we neglect magnetic polarization, which can be 
important for large, highly conductive particles.^ 



III. METALLIC PARTICLES 

This case is representative for absorbing particles. At 
low photon frequencies u> below the interband transitions 
region, metals can be well described by the Drude model, 
characterized by a DC electric conductivity Go and a di- 
electric function e = 1 + i47T(7o/w!^For a spherical par- 
ticle of radius a, we have a w a 3 (e — l)/(e + 2), and 
consequently 



Im{a(w)} 3wa 3 /Anao- 



(6) 



For sufficiently small particles, absorption dominates 
over radiative corrections, so that we can overlook terms 
proportional to \a\ 2 in Eqs. (|3|-([5]). Then, we find the 
closed-form expressions 



prad = ha3 
° 607T 2 C 3 (7o 
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where the subscript D refers to the Drude model and 

6j = 2Trk B T j /h. 

Equations Q and |8]) show that vacuum friction is al- 
ways producing stopping (MQ < 0), whereas the balance 
of radiation exchange between particle and free space 
can change sign depending on their relative tempera- 
tures. The general trend of these expressions is shown 
in Fig. [2^b). At low O, the torque scales as f2, whereas 
a steeper il 5 dependence is observed at faster velocities. 
Interestingly, a nonzero torque M cc O 5 is predicted at 
Tq = T\ = 0, despite the axial symmetry of the particle. 



IV. EQUILIBRIUM TEMPERATURE 

The power absorbed by the particle in the form of ther- 
mal heating p abs can be obtained from energy conserva- 
tion, expressed by the identity —MQ = P rad + p abs , 
where the left-hand side represents mechanical energy 
dissipation (stopping power). Using Eqs. §f§ and QSj> , 
we find 
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FIG. 2: Equilibrium temperature and stopping of a metallic 
sphere, (a) Normalized particle temperature at equilibrium 
(Ti/To) as a function of Q/6o, where 9o = 27rfc_gTb/?i (see Fig. 
[l]). (b) Universal normalized stopping power both at equilib- 
rium temperatures (solid curve) and at equal temperatures 
(To = Ti, broken curve). 



The particle equilibrium temperature is determined by 
the condition p abs = 0, and it is stable because 
dP ahs /dTi < [this inequality is obvious from Eq. 
but it can be easily derived in the general case from Eqs. 
([3]) -dSt]. Unlike conventional friction of a spinning object 
immersed in a fluid, vacuum friction is not always lead- 
ing to particle heating, as shown in Fig. [2^a) from the 
solution of Pg bs = 0. Actually, T x < T for finite tem- 
peratures and rotation velocities below 51 — 8q, whereas 
particle heating occurs at higher Q. The crossing point 
between these two types of behavior is independent of 
particle size a and conductivity ctq- 



At T = 0, we find 6 1 « 0.867 0, so that the ft 5 
dependence of Mu is maintained with the particle at 
equilibrium temperature. The loss of mechanical energy 
is then fully converted into a radiated power P™ d w 
0.013 ha 3 n 6 /c 3 a - 

It should be noted that having the particle at equi- 
librium temperature or at the same temperature as the 
vacuum results in significant differences in the stopping 
power [Fig.^b), calculated from Eqs. ([H])-©]. 
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FIG. 3: Power spectrum dP™ d /doj = ftoj[F(cj) - F(-lj)] [see 
Eq. Q] radiated by a metallic spinning particle for various 
rotation frequencies. Solid curves: emission at equilibrium 
temperatures. Dashed curve: absorption by a particle at rest 
and Ti = 0. The emitted-photon frequency u) is normalized to 
6*o = 2irkBTo/h. The inset shows the frequency of maximum 
emission at equilibrium as a function of Q/60. 



V. EMISSION SPECTRA 

The probability of emitting photons at frequency uj is 
given by T(u) — T(—lj) [see Eq. Q], which is normalized 
per unit of emission-frequency range. The emission pro- 
file at low rotation velocities (fl = 0.05 9q curve in Fig. [3]) 
mimics the absorption spectrum from a static cold par- 
ticle (dashed curve), also peaked around fkj ~ bksTi for 
Drude spheres. However, the maximum of emission is 
driven by Q for faster rotations (see inset and f2 = 5 9q 
curve in Fig. 3]), thus signalling a significant departure 
from standard black-body theory. 



VI. STOPPING TIME 

At low rotation velocity and finite temperature, the 
frictional torque acting on a metallic particle is propor- 
tional to fi [see Eq. ||8J)]. The correction to the particle 
equilibrium temperature [6\ s» #0 — (7/15)£l 2 /0 o ] can be 
then neglected to first order in f2, so the torque becomes 
M « -/3fi, where (3 = fta 3 6>o/307r 2 c 3 cr - From Newton's 
second law, we find an f2(t) = S7(0) exp(— i/r) time de- 
pendence of the rotation velocity, where r = I//3 is the 
characteristic stopping time and I is the moment of in- 
ertia. For a spherical Drude particle, we find 



(?ic) 3 - 2 



pa ctq 
(k B T )*- 



(10) 



where p is the particle density. 

Graphite particles are abundant in interstellar dust pi 
so we focus on them as an important case to study the 
rotation stopping time. The frequency-dependent dielec- 
tric function of graphite is taken from optical dataf^ 
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FIG. 4: Characteristic stopping time of spinning graphite par- 
ticles as a function of vacuum temperature. Solid curves: 
full calculation using measured dielectric functions for the 
graphite particles. 1 '' Broken curves: analytical Drude approx- 
imation [Eq. ( 10 1] . Various particle sizes and shapes are con- 
sidered: spheres of radius 10 nm and 100 nm, and an oblate 
ellipsoid of radius 10 nm and aspect ratio ij = 0.2. Low rota- 
tion velocities SI <C fesTo/?i are assumed (e.g., Q <C 21 GHz 
at To = IK). 



which could be eventually examined through measure- 
ments of rotational frequency shifts!^ Besides, relatively 
small stopping times are predicted for graphite nanopar- 
ticles, which ask for experimental corroboration (for ex- 
ample, using in- vacuo optical trapping setups). By anal- 
ogy to the Purcell effect the frictional torque can be 
altered due to the presence of physical boundaries that 
modify the density of states appearing in Eq. Q, thus 
opening new possibilities for controlling the degree of fric- 
tion (e.g., the torque can be strongly reduced at low tem- 
perature and small rotation frequency by placing the par- 
ticle inside a metallic cavity, which produces a threshold 
of p° in uj). 
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Appendix A: The fluctuation-dissipation theorem 



tabulated for different particle sizes, which differ due 
to nonlocal corrections. The \ow-lo behavior is well 
approximated by the Drude model with <7o = 2.3 x 
10 4 (2.0 x 10 5 )f2 _1 m _1 for spherical particles of radius 
a = 10(100)nm, where the response has been averaged 
over different crystal orientations. Plugging this into Eq. 
(10), we obtain the results shown in Fig. [4] by broken 
lines. Interband transitions become important in the re- 
sponse of graphite at frequencies above Tllo ~ 10 -2 eV, 
so we expect a deviation from Drude behavior at tem- 
peratures above ~ 100 K in this material. This is indeed 
confirmed by numerically integrating Eq. ^ with the full 
tabulated response of graphite to obtain r (Fig. |4j solid 
curves). For the particle sizes under consideration, stop- 
ping times are small on cosmic scales within the plotted 
range of temperatures, which are often encountered in 
hot dust regions.^ In cooler areas (To = 2.7 K), 100 nm 
graphite particles have a stopping time r ~ 0.6 billion 
years. 

Dust particles can adopt non-spherical shapes. In par- 
ticular, for oblate ellipsoids Eq. ^ [Im{a(uj)}} must be 
corrected by a factor ?y/9L 2 , where r\ is the aspect ra- 
tio (see inset in Fig. [4j and L is the depolarization fac- 
tor for equatorial polarization, approximately linear in 
^ED Also, / is linear in 77, thus leading to a r oc rj 2 de- 
pendence for fixed radius. We show in Fig. [4] the case 
a = 10 nm and rj = 0.2, which exhibits a significant re- 
duction in r compared to spherical particles of the same 
radius. In a related context, translational motion leads 
to thermal dragf^only at nonzero temperature and with 
similar stopping times. 

VII. CONCLUDING REMARKS 



The fluctuation-dissipation theorem (FDT) is a well- 
established result, first formulated by NyquispSl and later 
proved by Callen and Welton. - 1 it connects the fluctua- 
tions of the product of two operators with the dissipation 
expressed through the imaginary part of their response 
function. We give a simple derivation here, with a no- 
tation that is appropriate to deal with fluctuations of 
vacuum fields and particle polarizations. 

Let us consider a Hamiltonian perturbed by a term 

H'(t) = -<p(t)q(jt), 

where tp(t) is a time-dependent function, q(t) is 
an operator in the Heisenberg picture, related to 
its Schrodinger representation qg through q(t) = 
exp(iHot/h)qs exp(— iHot/h), and H$ is the unper- 
turbed Hamiltonian. In the Heisenberg representation, 
Schrodinger's equation becomes H'\<j>) — ihd\(f>)/dt, and 
we have \(f>) = ex-p(iH t /h)\<fis) . 

Under the condition H'(t) t ^ 0, the eigenstates of 
the total Hamiltonian become 



\<f> m {t)) = |m)- 



difH'{if)\<j> m {t')) 



dt'H'(t')\m), 



where the last line corresponds to first-order perturbation 
theory, and |m) is a state of the unperturbed Hamiltonian 
with energy E m (i.e., H \m) = E m \m)). 



The present results can be relevant to study the dis- 
tribution of rotation velocities of cosmic nanoparticles, 



The expected value of another operator p(t) is simply given by 



E< 



-E m /k B T 



[m\p(t)\m) + - I dt'<p(t')(m\[p{t),q(t')}\m) 



(Al) 



where 



Z = J2< 



-E m /k B T 



is the partition function at temperature T. The first term in Eq. (All reduces to (m\p(t)\m) = (m\ps\m), and from 
here, we can recast (|Alj) as 



(5p(t)) = (p(t) - #-00)) = / dt' X (t - t'Mt') 



where 



X (t - t>) = U{t f)i E e ~ Bm/fcBT H0(t), g(f)]|m) 



m 

is a susceptibility function. Now, using the closure relation |n)(n| = /, we can write 



(A2) 



\p(t),$(t')}\m)=J2[{m\ps 



n) {n\q s \m) e 



i(E m -E n )(t-t')/h 



(m\qs\n) (n\p s \m) e 



-i{E m -E n ){t-t')/n 



Inserting this back into Eq. ( A2 ), taking the time Fourier transform to work in frequency space, and using the identity 



we find 



X (w) = / dt X (t) e wt 
-1 



^2 (m\p s \n) (n\q s \ 



e ~E m /k B T _ e -E n /k B T 

hu! + E m - E n + iO+ 



Incidentally, the zero-temperature susceptibility reads 

(0\ Ps \m) (m\q s \0) (0\qs\m) (m\ps\0) 



T-vO 



E 



_ Tiuj + E - E m + iO+ Huj + E m - E + iO^ 
Finally, the dissipation associated to \ can be written 

Im{ x H}=(l- e ^M ZLj%-iW*BT \ m \p s \n) {n\q s \m) 8{hw + E m -E n ). 



(A3) 



Similarly, we can write the average over fluctuations as 

,-E m /k B T j(E m -E n ){t-t')/n 



S(t-t')^(p(t)q{t')) = \Y. e ~ : ' ' h <>»\i*\»)<>'\<is\>'>) 



duj 
2n 



■j(t-t') 



where 



2?r?i 



S(u>) = J2^ Em/kBT HPs\n) (n\q S \m) S(hu + E m -E n ). 



(A4) 
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The relation between S(ui) and Im {x(w)} that one obtains by comparing Eqs. (A3 ) and ( A4 1 constitutes the general 
form of the fluctuation-dissipation theorem: 



S(u) = 2% [n(u) + 1] Im{x(w)} , 



(A5) 



where 



n(oj) 



1 



e huj/k B T _ I 



is the Bose-Einstein distribution function. 

We can formulate a more useful relation by noticing that (p(t)q(t')) is a function of t — t' , so that its double Fourier 
transform satisfies 

(0(w)$(w'))= / dtdt' e lLOt+lu,t 'S{t-t')= [ dre lUT S(T) f dt' e^ +LU ' }t ' = 2ir 5{lo + J) S(u), 



(A6) 



and from here we find the expression 

(p(w)g(a/)) = 4nh [n(w) + 1] Im{x(w)} 5{uj + u)'). 



Proceeding like abov e, th e Fourier transform of the fluctuation (q(t')p(t)) reads (q(uj')p(uj)) = exp(— Tu^/feeT) S(w), 
which together with ( A5 ) leads to 



(q(ui')p(uj)) = Anh n(uS) Im{x(w)} 8{u> + a/). 



(A7) 



Finally, it should be noticed that p{uj)q(ui') is not an observable in general, but the symmetrized product is Hermitian, 
and therefore, an observable. From Eqs. (A6| and (A7|, we find 



1 



(p(w)g(w') + q(u')p(u)) = inh 



n(oj) + 



Im {x(w)} S (u + L)'). 



(A8) 



Equations (A6)-(A8) are general forms of the FDT. Next, we formulate specific applications for dipole and electric-field 
fluctuations. 



1. FDT for fluctuations of the dipole moment 



Now, we can apply the above general expressions of the FDT to dipole-dipole fluctuations, with the identifications 



p(t) -> Pi (t), 
q{t) Pj(t), 
x (t) -> atij(t), 
<p(t) ^ 

where pi and pj are components of the dipole moment along directions i and j, respectively, Ej is the electric field 
along j at the position of the dipole, and is the component of the polarizability tensor. The interaction 

Hamiltonian is H'(t) = —Ej(t)pj(t), where Ej is regarded as a time-dependent function and qj as an operator. The 
susceptibility acts in frequency space according to (6pi(uj)} = aij(u>)Ej(uj). With these substitutions, the FDT [Eqs. 
(A6)-(fdt3)] takes the forms 



(Pi(w)Pj( w ')) 
(Pj(uj')pi(u})) 

(pi(u})pj(u)') +pj(u)')pi(u))) = AitTi 



4irh [n(w) + l] Im{o!y(a;)} 5(u}+u/), 
Airh n{uj) Im {otij(w)} 6{u + u), 

n(uj) + - lm{aij(ui)} 8(ui + ut'). 



(A9) 
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2. FDT for fluctuations of the electric field 



Similarly, the fluctuations of the electric held can be analyzed with the substitutions 



p(t) -> Ei{v,t), 
q(t) ->■ E,{r',t), 

where J3j(r, t) and Ej(r',t) are components of the electric held along directions i and j at positions r and r', respec- 
tively. The Green tensor of the electromagnetic held in vacuum is given, in frequency space w, by 



Gy(r,r» = 



exp(ikR) 



(k 2 R 2 + ikR - 1) Sij - {k l R 2 + 3ikR, - 3) 



RiRj 



(A10) 



where R = r — r' and k = uj/c. The interaction Hamiltonian is again H' — —pjEj, but now pj is a time-dependent 
function and Ej is an operator. The susceptibility acts as 

(5E i (cj))=G ij (r,r',uj)p j (uj), 

which is consistent with the definition of Eq. (A10). 2 ^The FDT takes the forms 

(E i (r,u))E j (r',Lj')) = iirh [n(w) + 1] ImjGy^r', w} 5(u + u'), 
(E j (i',u')E i (T,u)) = 47rhn(u;)lm{G tJ (r,r',u;)} 8{u+u/), 

1 



(E l (r,u J )E J (r',uj') + E :i (r',uj')E t (r,uj)) = AttH 



n(w) 



Im{G y (r,r',tj)} S(lj + lj'). 



Interestingly, for r = r' we have Im{Gij(r,r, u>)} = (2ir 2 ui/3)p°Sij, where p° is the photonic local density of states 
(p° = cu 2 /ir 2 c 3 in vacuum), which yields 



{Ei[T,u)Ej[T,J) +E J (r,w , )E l (r,oj)) = 



8ir 3 hojp° 



n(tj) + 



6^ 5(u + (J). 



(All) 



Appendix B: Polarizability of a rotating particle 

Rotational motion affects the polarizability a of a small spinning particle. This quantity represents bubble diagrams 
made up of virtual (de-)excitations. The corresponding matrix elements involve particle ground and excited states of 
electronic and vibrational nature, which do rotate with the particle. Therefore, a provides the relation between fields 
and dipoles expressed in the rotating frame. The external electric held has to be transformed to the rotating frame 
before multiplying by a to obtain the induced dipole, and this in turn has to be transformed back to the rest frame. 

We assume rotation around the z axis, and denote by x and y the unit vectors in the lab frame, and x' and y' the 
unit vectors in the rotating frame. With the particle rotating with angular frequency the angle between x' and x 
increases with time as fit. Then, expressing the external electric held in frequency space w, we can write 

E(t) = J ^'^[E^i + Ey^y] 

= J ^e-«K( W )x' + £;Hy'], 

The relation between the held in the lab and rotating frames can be obtained by writing x and y in terms of x' and 
y', and then absorbing the fit dependence in the time exponentials within the above integral. We obtain 

E' x {lj) = ~[E x (w + ) + E x (w-)-iEy(u + )+iEy(u-)], 
E' y {u) = ^[iE x (u> + )-iE x (Lj_)+E y (Lj + )+E y (u_)}, 
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where w± = w±fl. Besides, the induced dipole in the rotating frame is p' x (oj) = a_\_{uS)E' x {u>) and p' y (uj) = a±(u))EL(u), 
where a± is the static-particle polarizability for polarization perpendicular to the rotation axis. We assume the particle 
to have axial symmetry around z, so that polarization along that direction, described by ay, is unaffected by the 
rotational motion. Finally, we need to transform the induced dipole from the rotating frame to the lab frame. 
Proceeding like above, we have 



P V (UJ) 



[Px(u+) +P x {u-) + ip' y (u+) - ip'yi^-)] , 
[-ip x (u> + ) + ip' x (w-) +p' y (u> + ) +p' y {uj-)\ 



Putting all this together, we find the induced dipole in the rotating frame to have the same frequency as the external 
field for this particular case of axially-symmetric particles. More precisely, the relation between electric field and 
induced dipole in the lab frame is described by an effective polarizability a eS (uj) according to 







P y(uj) 




Pz{u)_ 







<( 


















E x {u) 







E y (u) 






E z (u>) 



where 



ffiu) = -[a ± (u, + tt)+a ± (u)-!l)], 



■[a±(u + D)-a±(u-n)], 



(Bla) 

(Bib) 
(Blc) 



Like a, this effective polarizability satisfies the retarded response condition a cff (— uS) = [a cff (w)]*. 



Appendix C: Derivation of the vacuum torque 

The torque exerted by an electric field E on a dipole p is given by 

M = p x E. 



(CI) 



This is intuitively understood by considering the dipole to be formed by a negative point charge located at the origin 
and a positive point charge on which the electric field produces a force; the position vector of the positive charge 
enters both the dipole and the torque, thus yielding Eq. (CI I. The vacuum torque can be obtained from (CI), with 



p and E produced by polarization and field fluctuations, respectively. Noticing that these two kinds of fluctuations 
are uncorrelated, the vacuum torque for a particle with symmetry of revolution around its rotation axis z is given by 



71/ 



z • (p fl x E 

/rfl Find „fl pii 
{Px^y ~Py h x 



ind 



a x E a) 

Px d E S y - 



El), 



(C2) 



where the average is taken over fluctuations of the vacuum electric field E fl and the particle polarizability p fl . 

With the particle at the origin, we obtain the induced field in frequency space from E md (r, uS) — G(0, r, lj) ■ p (w), 
where the component s of G are given by Eq. ( A10 ) and the r — > limit is to be understood in what follows. Similarly, 



the induced dipole in (C2| is obtained from p im (cj) = a cS (ui) -E fl (0, ui). Writing the fields and the dipoles in frequency 



space and using these equations for the induced fields and dipoles, we obtain from ( C2 ) 
dujduj 



M = 



G yx (0, r, W ) (pI{uj)pI{u')) + G yy (0, r, w) 

-G w (0,r, W ) (pg(w)pg(a/)) - G I9 (0,r, w ) (pg(w)^(a/)) 
+af y (io) (£"(0,w)£ 9 9 (0,w')) + af(w) (^(0, w)^(0, w')) 

(w) (E*(p,u>)E%(0, U ')) - af y (u,) (^(0,w)^(0,a/)> 



v cir 

i/.r 
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Explicit expressions for G and a cB are taken from Eqs. (A10) and (Bl I, respectively. Finally, we apply the FDT from 
Eqs. ( A9 ) and ( Al 1 1 , assuming that the products of the p and E observables are symmetrized as explained at the end 



of the FDT section. The field fluctuations occur in the vacuum at a temperature T , so that the corresponding FDT 
yields a Bose-Einstcin distribution at that temperature, which we denote hq. Likewise, the polarization fluctuations 
take place in the particle at temperature T±, and thus one obtains a distribution ri\. After some straightforward 
algebra, we find that only the imaginary part of G survives in the above integral. Finally, the torque reduces to 



M 



where 



r(w) = (2nojp / 3) 2Im{a ~ fl)} [m(w - 0) - n {uj)] 



(C3) 



(C4) 



is a spectral distribution function (see below). A similar result is obtained for particles without axial symmetry, but 
with orthogonal axes of polarization, and rotating around one of them, z; in this case, after time averaging, we find 
that a± has to be substituted by the average of the polarizability over directions perpendicular to the rotation axis, 



*;/;/ 



)/2- 



A similar derivation can be carried out for the torque produced by magnetic polarization, which is negligible for very 
small particles in which the magnetic polarizability is generally small, although it can play a role for metallic particles 
of large conductivity, leading to imaginary parts of the electric and magnetic polarizabilities that are comparable in 
magnitude (the magnetic effect can be even dominant). 



Appendix D: Derivation of the radiated power 

Following the intuitive explanation of the torque acting on a dipole (see previous section) , we can argue that the 
work exerted on a dipole by an electric field is P = E • dp/dt. Here, we are interested in the radiation produced by 
friction acting on the particle, or equivalently, the work done by the particle on the vacuum, which we can write in 
terms of fluctuating and induced dipoles and fields as 

prad _ _ ^ E i„d . Q p fl/ dt + E fl . dp^/Qt) . 

After lengthy but straightforward algebra, mimicking the steps followed in the calculation of the vacuum torque, we 
obtain 



/oo 
-oo 



where T(ui) is the same as in Eq. (C4). Now, we can interpret V as the spectral distribution of the rate of emission 
(when ojT > 0) or absorption (ojT < 0). It should be noted that u) in these expressions has been carefully ensured to 
correspond to the photon frequency in the lab frame. 

Incidentally, we are missing in T the contribution of polarization along the rotation axis z, which can be easily 
derived from half the contribution of radial polarization for a non-rotating particle. The full result becomes 

r(w) = (2trV73) |2Im{a±(o;-n)} [m(a; - fi) -n Q {u))] +Im{a||(w)} [m(w) - no(w)] }. (Dl) 

Actually, the new term inside the curly brackets, together with the leading prefactor, yield an odd function of w, so 
that this term vanishes under the integral of Eq. (C3), and therefore, we can use Eq. (Dl) both for the torque and 
for the radiated power. 



Appendix E: Quantum-mechanical coupling of rotational motion and free photons 

In a quantum-mechanical description of rotational vacuum friction, the interaction between the particle and the 
free-space electromagnetic field can be described by the HamiltoniarP^ 

H'=iJ2\f*^(<l-ei)l"j-*t]> ( E1 ) 

3 
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where the sum runs over photon states j of polarization vectors £j and frequencies Uj, V is the quantization volume, 
and dj are photon creation and annihilation operators, respectively, and d is the particle-dipole operator. We 
have considered the particle to be small compared to the wavelength of the involved photons, so that the use of the 
dipole approximation is justified. 

The rotational state of a spinning particle can be described as a combination of eigenfunctions exp(imtp) / \/2n, 
where tp is the rotation angle and m is the azimuthal quantum number. For rotation velocity f2, the values of m are 
peaked around m ~ I£l/h, where I is the moment of inertia. The angle tp enters Eq. (El ) through the transformation 



of the dipole operator from the lab frame (d) to the rotating frame (d'). More precisely, 

d x = d' x cos tp — d' sin tp, 



d y = d' x sirup + d' cos tp. 



The terms in sincp and costp can produce transitions m — > m ± 1, while d' x and d' y generate internal transitions in 
the particle. This is accompanied by the emission or absorption of one photon, according to Eq. (El). In other 



words, rotational friction is associated to transitions involving simultaneous changes in (i) the rotational state of the 
particle, (ii) its internal state, and (iii) the number of photon states. All of these elements must be present in the 
interaction described by Eq. (El). This implies that non-absorbing particles (at least within the relevant range of 
photon frequencies) cannot undergo friction, because they do not possess the excited states that are necessary to 
sustain internal excitations. 
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